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Abstract 

The Kadomtsev-Petviashvili I (KPI) is considered as a useful laboratory for 
experimenting new theoretical tools able to handle the specific features of integrable 
models in 2 + 1 dimensions. The linearized version of the KPI equation is first 
considered by solving the initial value problem for different classes of initial data. 
Properties of the solutions in different cases are analyzed in details. The obtained 
results are used as a guideline for studying the properties of the solution u{t, x, y) 
of the Kadomtsev-Petviashvili I (KPI) equation with given initial data u(0, x, y) 
belonging to the Schwartz space. The spectral theory associated to KPI is studied 
in the space of the Fourier transform of the solutions. The variables p — {pi,P2} 
of the Fourier space are shown to be the most convenient spectral variables to 
use for spectral data. Spectral data are shown to decay rapidly at large p but to 
be discontinuous at p = 0. Direct and inverse problems are solved with special 
attention to the behaviour of all the quantities involved in the neighborhood of 
t — and p = 0. It is shown in particular that the solution u{t, x, y) has a time 
derivative discontinuous at t = and that at any t ^ Q \t does not belong to the 
Schwartz space no matter how small in norm and rapidly decaying at large distances 
the initial data are chosen. 
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1 Introduction 



We consider the initial value problem for the Kadomtsev-Petviashvili equation ||, |j 
in its version called KPI 

(ut - 6uux + u^xx)x = iuyy , u = u{t, X, y) , 

■u{0,x,y) = u{x,y) , (1.1) 

for u{t,x,y) real. Already in 1974 Q it has been acknowledged to be integrable since 
it can be associated to a linear spectral problem and, precisely, to the non-stationary 
Schrodinger equation 

^_,d^ + dl-u{x,y))'P = 0. (1.2) 
However, to building a complete and coherent theory for the spectral transform of the 



potential u{x,y) in (1.2) that could be used to linearize the initial value problem of (1.1) 
resulted to be unexpectedly difficult. The real breakthrough has been the discover that 
the problem is solvable via a non local Riemann-Hilbert formulation H, ||] . 

Successively other progresses have been made. The characterization problem for the 
spectral data was solved in . The extension of the spectral transform to the case of a 
potential u(x,y) approaching to zero in every direction except a finite number has been 
faced in [||. The questions of the associated conditions (often called 'constraints') and 
how to choose properly in the evolution form of KPI 

ut{t, X, y) - 6u(t, X, y)uxit, x, y) + u^xxit, y) = 39^^ Uyy{t, x, y) (1.3) 

have been studied in H. Some relevant points on constraints were known earlier, but 



have not been published [Q. See also ||TT|. In the article |12 some rigorous results on 
the existence and uniqueness of solutions of the direct and inverse problem have been 
proved. 

In this article, as well as in the previous one ||l^ by the same authors, the description 
of the KPI equation by means of the inverse spectral transform is revisited. We think 
that this is necessary since some specific properties of the solutions of the initial value 



problem (IT) and, more generally, of the 2 + 1-integrable equations are not properly (or 
in any case not exhaustively) described in the literature. Moreover, these properties are 
specially interesting since they seem to be peculiar of the integrable equations in 2 + 1 
dimensions. Our study is essentially based on some new theoretical tools in the theory 
of the spectral transform developed in ^, [l^, [l5|| . 



In this article we consider the equation (1.1) for initial data u{x,y) belonging to the 



Schwartz space S. The main results we get as regards the properties of the solutions 
and the theory of the spectral transform can be summarized in the following points (the 
integrations when it is not differently indicated are performed all along the real axis from 
—CO to +co). 

1. The solution u{t,x,y) immediately leaves the Schwartz space since at any time 
t ^ develops a tail slowly decreasing for x too according to the following 
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formula 

u{t,x,y)^ I ldx'dy'u{x\y')+o{\x\-^'^). (1.4) 



A.Tl\/?)tx\x\ 

Only in the opposite direction, i.e. for x — > —too, it decreases rapidly. 



2. The problem of the proper choice of d^^ in (l.S) requires a special investigation in 
the neighborhood of the initial time i = 0. The final answer is that, for initial data 
u{x, y) £ S which are arbitrarily chosen and, therefore, not necessarily subjected 
to the constraint 

dx u{x, y) — , (1-5) 



the function u{t,x,y) reconstructed by solving the inverse spectral problem for 
( L2 ) evolves in time according to the equation 

ut{t, X, y) - 6u{t, x, y)ux{t, x, y) + u^xxit, x,y) = 3 dx' Uyy{t, x' , y) . (1.6) 



Notice that the possibility that the inversion of dx can depend on t was already 
mentioned in the literature for the Davey-Stewartson equation pGf . 

3. The time derivative Ut{t,x,y) has at t = different left and right limits and the 
condition 

dx u{t, x,y) — (I-''') 



is dynamically generated by the evolution equation at times t ^ 0, i.e. for these 
times we recover the result obtained in ||^. 

4. The primitives or their antisymmetric combination ^ ^ J-oo ^ cannot 

be exchanged with the limit t 0. In fact we have 

/X f + OG 

dx' u{t,x'y) = — lim / dx' u(t, x' ,y) — 
-oo *-+0i. 



1 



+ 00 



and 



4^1^02'' ~' )^2:'M(t,a;',y) = / dx'u{x',y) 



/X r + 00 

dx' u{t, x' ,y) = — lim / dx'u{t,x',y) 
-00 ^^^^Jx 



tll™o2W ~J Wu(i,a;',2/)--y dx' u{x' ,y) . (1.9) 

Thus the two forms of writing down the evolution form of the Kadomtsev-Petviashvili 
equation suggested in and in (1.6) are equivalent but only in the last formulation 
the integration over x and the limit t — s- can be exchanged showing explicitly the 
discontinuity at t = 0. 
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5. Conditions ('constraints') higher than (1.7) can be obtained. But some special 
care is needed due to the asymptotic behavior of u as given in (|l.4|). E.g. in the 
condition 

dxxuy{t,x,y) = , t^O, (1-10) 



/ 



the y-derivative cannot be extracted from the integral. Higher conditions arc non- 
hnear and can be given by means of recursion relations. 



6. In spite of the dynamic condition (1.7) we have that at all times 

dx dyu{t,x,y) = / dxdy u{0,x,y). (l-H) 



The integrations in the l.h.s. are performed in the order explicitly indicated. Since 
no any constraint is imposed to the initial data this quantity is not necessarily 
equal to zero and gives a nontrivial integral of motion. 



7. The properties of the Jost solutions of the equation (1.2) essentially differs at t — 
and t ^ 0. In particular, at t = the coefficients of the asymptotic 1/k-expansion, 
where k denoted the spectral parameter, depend on the sign of the imaginary part 
kcj (cf. ijl^ and [|l5j), while at t ^ these coefficients in the limit i — > depend on 
the sign of t. 

8. Spectral data do not belong to the Schwartz space S and no any small norm or 
smoothness condition can improve this situation. 

All along the article we find convenient to consider the problem in terms of the Fourier 
transform of u 

v{t,p) = v{t,p,,p2)^j^ JJdxdye'^P^--P^y^u{t,x,y), (1.12) 

where p denotes a 2-component vector 

p={pi,P2}. (1.13) 



Then, the initial value problem (1.1) for the KPl equation can be rewritten as {d^q = 
dqidq2) 

pi — — = -3ipi jd^qv{t,q)v{t,p- q) ~ i {pi + 8^2) v{t,p) , 



dt 

viO,p)^vip)eS, (1.14) 



where v{p) is the Fourie r tr ansform ( 1.12 ) of the initial data u{x,y). The problem of 
properly defining d^^ in (^]^) is reformulated as the problem of properly regula rizi ng the 
distribution Pi^. In terms of the Fourier transform v(t,p) we have instead of (|l.6| ) 

dv{t,p) ^ [d\v{t, q) v{t,p- g) - » ^ v{t,p) (1.15) 



dt -^..y^ p^ + m 
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and condition (^]^) takes the form 



lim i;(i,pi,p2) = 0, for t ^ 0, (1.16) 

where the limit is understood in the sense of distributions in the p2 variable. 

According to the usual scheme of the spectral transform theory we expe ct t hat the 



spectral data evolve in time as the Fourier transform of the linear part of This 
makes clear that the special behavior in time of the solutions of KPI is just complicated 
by the nonlinearity of the KPI but is, in fact, inherent to the singular character of its 
linear part. Due to this we start the article by considering in section 2 the linearized 
version of the KPI equation. The above mentioned properties are proved for the class of 
Schwartz initial data. As well some more general classes of initial data are considered, 
for which the condition of unique solvability of the initial value problem is preserved 
and the corresponding evolution forms of the equation are presented. The study of the 
nonlinear case for these initial data is an open problem. It is show that in spite of the 
discontinuity at < = the linearized evolution equation for initial data belonging to the 
Schwartz space is an Hamiltonian system. Thus it is natural to expect that the KPI 
equation with this type of initial data is Hamiltonian too. But since a formulation of the 
simplectic and Poisson structures in terms of spectral data meets with some difficulties 
even in 1 + 1 dimensions, i.e. in the KdV case (see, e.g. (l^), we postponed this problem 
to a forthcoming publication. 

In section 3 the inverse spectral transform method for the KPI equation is revisited. 
We use a specific form of the equations of the inverse problem suggested in The 
advantages of this formulation as well as the connections with the results obtained in 
are given in section 4. A lemma and some propositions on distributions that are of 
common use in the article are reported in the appendix. 



2 The Linearized Equation 

We consider the linear part of the KPI equation 

dAdtUit, X, y) + dlU{t, X, y)) = idlU{t, x, y). (2.1) 

and its initial value problem, i.e. we are searching special classes of solutions U{t,x,y) 
satisfying 

UiO,x,y) = Uix,y). (2.2) 

The solution U{t,x,y) and the initial value U{x,y) can be distributions in the (t,x,y) 
variables belonging to the space S' dual to the Schwartz space S, but U{t, x, y) must be 
continuous at t = in order to have a meaningful initial value problem. 

Since ( |2.l| ) is not an evolution equation its initial value problem is singular and 
we must pay special attention to the behaviour of the solution U{t,x,y) at the initial 
time t = Q. From the study of this problem we expect to get useful guidelines for 
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the KPI equation. In fact, the additional nonhnear term present in the KPI equation 
comphcates the analysis but does not change the specific behaviour of the solution in the 
neighborhood to the initial time. 



Performing the Fourier transform of the equation (2.1) we get 

dVit,p) 



Pi 



dt 



= -i{p\ + 'dpi)V{t,p) 



where 



V{t,p) = V{t,p^,P2) 



1 



(2^)2 



dxdye'''P'''-P'-y^U{t,x,y) 



and p denotes a 2-component vector 

P= {Pl,P2}- 

For any Pi ^ the general solution of ( ^.3[ ) satisfying the initial value requirement 

ViO,p) = V{p) 

is given by 

V{t,p) = exp ( —it 

We choose V{p) to be the Fourier transform of the initial data U{x,y) 

1 



Pi 



V{p), 



Vip) 



i2n) 



dx dy e 



ipix-ip2y 



U{x,y). 



(2.3) 
(2.4) 

(2.5) 

(2.6) 
(2.7) 

(2.8) 



Then (2.6) is the Fourier transform of (2^) and the initial value problem defined in (2.1) 
and (2.2) has been transformed into an equivalent one in the Fourier transformed space 
of variables p. 

A natural regularization at pi = of the distribution in (2.7) is just formula (2.7) 
itself considered for all values of the variables (t, p) including pi ~ 0. In order to compute 
its time derivative it is convenient to consider it as the limit in the sense of S' 



Vit,p) 



lim exp 



-it 



Pi 



3pi 



Pi + iet 



V{p). 



(2.9) 



Then we derive with respect to t both sides and exchange the time derivative and the 
limit in the r.h.s., which is allowed in the space S' since a test function multiplied by 

exp ^— li^^:^^ still belongs to the Schwartz space. We get (see lemma A.l in the 

appendix for an alternative proof) that V(t,p) satisfies the evolution equation 



dV{t,p) ^ _.p\ 



■ipl 



dt 



pi + iOt 



V{t,p). 



(2.10) 
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Note that the special regularization of 1/pi obtained in the r.h.s. is fixed by the re- 
quirement that V{t,p) evolves in time according to (2/7) for all values of pi including 
Pi=0- 

The above considerations can be extended to the cases in which sgnpi V{p), P2V(j>) 
and belong to S. They need to be considered in details. For instance in the 

first case the time derivative and the limit cannot be exchanged any more and the time 
derivative of V{t,p) is not defined at i = 0. A complete study of all these cases is 
considered in the following. 

We conclude that V{t,p) furnishes via the Fourier transform 

Uit,x,y)^ jifpe-'P'-'+'P^yVit.p) (2.11) 

a solution [/(t, x, y) belonging to the space of distributions S' satisfying the initial data 
equation (|2.2| ) and the evolution equation 

dtU{t,x,y) + dlU{t,x,y)^3 f dx' d^yU{t,x' ,y). (2.12) 

J —too 

Note that this equation must be considered as an equation for distributions in the vari- 
ables X, y and t. In particular at the initial time i = it can be continuous or discon- 
tinuous with left and right limits or not defined at all. In the following sections we shall 
give an explicit example for all these three possibilities. 

Since the most general regularization of V(t,p) at pi = can be obtained by adding 
a distribution with support the point pi — (see [^) we conclude that the general 
solution of the problem (2.1), (2.2) can be written in the form 

U{t,x,y) + nit,x,y), (2.13) 

where 

r^(t,2:,y)= [d'^pe-'P^^+'P'yuj{t,p) (2.14) 



is the Fourier transform of a distribution uj{t,p) concentrated at pi = 0. We have 
therefore to find the general solution of the problem 

Pi^^^^ = -i{pt + 3pl)Luit,p), uj{0,p) = 0, (2.15) 
where uj{t,p) is a finite sum of the (5-function 5{pi) and its derivatives of the form 

N 

c.(t,p) = ^<5(")(pi)c.„(i,p2), (2.16) 



with ujn some distributions in the p2 and t variables. By inserting ( 2.16 ) into ( 2.15 ) and 
by using the known properties 

Pi<5(pi) = 0, pi<5(")(pi) = -n5("-i)(pi), n=l,2,... (2.17) 
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we get 

pluJNit,P2) = 0. (2.18) 

and, then, by recursion 

{pl)"+'LON-n{t,P2)=0 n=l,...,iV. (2.19) 



Using again property (2.17) for p2 we get that a;„(<,p2) is a sum of the (5-function (5(p2) 
and its derivatives up to the order 2{N — rt) + 1 with time dependent coefRcients. 

This resuh can be re-expressed in terms of the x and y variables by saying that 
J7(t, X, y) is a polynomial in x and y of the form 

N 2(Af-ri) + l 

nit,x,y)^J2 E «n™(t)x"y'". (2.20) 

n— 7n—0 

The coefficients anm{i) obey an undetermined linear system of ordinary differential 



equations, which is obtained by inserting U,{t,x,y) into (2.1). It is easy to check (we 
omit for shortness the corresponding calculations) that all the coefficients anm{t) for 
m = 2, . . ., 2(7V — n) + 1 can be expressed as linear combinations of a„o(0 o,ni{t) and 
their time derivatives. These functions are arbitrary up to the condition that f2(f, x, y) 
is identically zero at the initial time t — 0. 
We conclude that 

Proposition 2.1 For the initial data U{x,y) we are considering, since they are vanish- 
ing at large x (or y), if there exists a solution U{t,x,y) of the linearized KPI equation 
vanishing at large x (or y) it is uniq ue. This solution is given by the Fourier transf orm 



of the distribution V{t,p) defined in (2/1) and satisfies the evolution equation { 2.1i ) 



In the following we will consider four different classes of initial data vanishing at large 



x 01 y and we will show that in three cases the Fourier transform of (2.7) satisfies (in 
the sense of the distributions) ( ^.12 ) at all times including t — 0, vanishes at large x or 



y and is therefore unique, while in one case the time derivative of the Fourier transform 
of ( ^.7| ) at f = is not defined and therefore the evolution form of the linearized KP 
equation has a (unique) solution in the space considered for any t ^ 0. 

2.1 Some special classes of solutions of the linearized equation 

We want to build special classes of solutions of the evolution version of the KPI equation 



in (2.12). This can be done by noting that, for any distribution G{t,x,y) satisfying the 



partial differential equation 

{dtd^-3dl)G{t,x,y) = 0, (2.21) 

the convolution 

T{t,x,y) = J ldx'dy'Git,x - x',y -y')Uo{t,x',y') (2.22) 
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of G with an arbitrary solution Uq of the equation 

{dt+dl)Uo{t,x,v)^0 (2.23) 
satisfies the linear partial differential equation 

{dt + dl) d^T = 3d^T (2.24) 



and that this last equation can be considered the time evolution version of (2^) for 

U{t,x,y) =d^T{t,x,y). (2.25) 
The general solution of ( 2.23| ) can be written as 

Uo{t,x,y) d^pe-'P'-+'P-y-''PWo{p) (2.26) 

with Vo{p) E S'. We require the following integrals 

fdxUo{t,x,y) =JdxUQ{0,x,y) = 2T:Jdp2e'P-yVo{0,p2) (2.27) 

fdyUo{t,x,y) = 27r/dpie-^f^^-^*P'Vo(pi,0) (2.28) 



to be well defined and convergent, i.e. we require Vo{p) to be at pi = and at p2 = a 
well defined distribution in the variable p2 and pi , respectively. 

Since we are studying an initial value problem special attention must be paid to the 
distribution G{t, x, y) and its derivatives at t = 0. In particular, since 

UiO,x,y)= f fdx'dy'd^Git,x-x',y-y%^oUoiO,x',y') (2.29) 



the distribution G(t, x, y) must be defined at i = 0. 
Let us, now, consider the quadratic form 

Q{t,x,y) = l2xt~y^ (2.30) 

and let us define 

= \Q\^^{-Q), Qi = \Q\^m)- (2.31) 

These distributions, following can be defined by analytical continuation in A. Both 
have simple poles at A = —1, —2, . . . , — fc, . . . and Q'^ has additional simple poles at A = 
—3/2, —3/2 — 1, ... , —3/2 — fc, . . . The residua at the first kind of poles are distributions 
with support on the cone surface Q = while the residua at the second kind of poles are 
distributions with support on the origin, i.e. on the vertex of the cone Q = 0. 

We are interested in using these di strib utions to build distributions G(t, x, y) satis- 
fying the partial differential equation ( 2.21 ). For all the distributions introduced in the 
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following we are using notations and definitions of jist . For instance the distribution ^ 
is defined by 

f I . f 1 

dx —ip(x) = \im / dx —Lp(x) 

J X ^ ' e^(>J\^\>^ X 

where ^p{x) is an arbitrary test function belonging to the Schwartz space S. 
We are considering four cases: 



(2.32) 



Case i) 



G+{t,x,y) 



(2.33) 



Case ii) 



Gl{t,x,y) 



2aS{x)i}{~ay), 



(2.34) 



Case iii) 



G-{t,x,y)^-Qy' 

TT 



(2.35) 



Case iv) 



G-{t,x,y) 



TT X 



(2.36) 



The corresponding solutions will be noted, respectively, U+,U1,U^,U-. The initial data 
U(0, X, y) can be computed in terms of the initial data J7o(0, a;, y) of the evolution equation 
( ^.23| ) that according to the requirements made for the integrals in (2.27) and ( ^.28 ) must 
vanish at large x and y. It results (see in the following) that ?7(0, a;, y) vanishes at large 
X and does not necessarily vanish at large y. This is sufficient to assure that all four 
cases furnish solutions belonging to the special class considered in the previous section 
that admits one and only one solution for the initial value problem. 

Case i) It is convenient to exploit the general theorem about tempered distributions 
assuming that any distribution can be expressed as derivative of a regular function. 
It is easy to verify that 



Q-^'^ - n5{y)d{xt) = I- 
dy 



arctan - 



y 



\Q\ 



1/2 



-^gnyU{Q) 



(2.37) 



and that 



7rd{y)d{xt) 



dy"^ 



y arctan ■ 



,1/2 



\y\]m) 



(2.38) 
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From (2.37) we have 



and from ( 2.3q ) by deriving with respect to x 



d_ 
dx 



Qj^l"^ - TrsgntSix)S{y) 



2 dy"^ [x 



1/2 



and, consequently, at t = 



— G+(t,x,y) 



= 5{x)S{y). 



t=±o 



From ( 2.40 ) by deriving with respect to t we get 



dt dx dy'^ 



27rS{t)S{x)S{y). 



From this equation, due to ( 2.41 ) , it fo llows that G+ satisfies ( [2.21 ) . 
We have therefore from ( 12.25^ and ( |2.4lD 



U+{0,x,y)^Uo{0,x,y) 



and from ( p^ ) and ( ps| ) 



lim r+(t,a;,y) = / da;'J7o(0, x', y). 



Case ii) For studying the distribution G"^ wc consider the identity 



^Q+^^\6nt6'iy)^(xt) = - 



From it at t = we derive 



9y2 



6t I arctan 7777 — — sgny ^{Q) 



X 



= 



t=0 



and, consequently. 



and 



Gl[t,x,y) 



t=o 



-2aS{x)'d{-ay) 



^G%{t,x,y) ^-2<j6'{x)d{-ay) 
^ t=o 



9x ^ 

By deriving ( ^.45 ) with respect to x we have 



d_ 

dx 



X 



-1/2 



1 92 



V ^1/2 



67r|i|<5'(y)5(:.) = 



11 



and, then, by multiplying by sgnt and by deriving with respect to t 
d d 



3 

dt dx dy^ 



and at t = 



[sgnt^Q;^/^] = -Q^5'{y)5{x). 
ataxi X i 

We deduce that G\ satisfies < ^.2l\j . 
We have therefore 

U^i0,x,y)^2r dy'dMO,^,y') 

J (TOO 

and y 

lim T^{t,x,y) = 2 [ dy'UoiO,x,y'). 

Note that the initial data satisfy the constraint 

fdxUl{0,x,y) = 0. 

— 1/2 

Case iii) Let us, now, consider Q_ . We have 



QZ 



1/2 



dy'^ 



1 \y-\Q\'^^\ 



and by deriving it with respect to x 



dx 2 dy"^ 



and at t = 



or 



ox 



t=0 



— G^{t,x,y) 



t=o 



= S{y). 

TTX 



(2.50) 
(2.51) 

(2.52) 
(2.53) 

(2.54) 



(2.55) 

(2.56) 

(2.57) 
(2.58) 



By deriving, in addit ion, with respect to t we obtain that G_ satisfies ( ^.21 ). However 
by considering ( 2.55| ) at i = it results that G_ is not defined at i = and, since G_ 
satisfies ( 2.21 ), that also dtdxG^ is not defined at t = 0. 
We have therefore 



U^O^x.y) 



1 fdx'^-^UoiO,x',y). 



TT I X — X 



(2.59) 



However, the limit T-(t,x,y) for t — > ±0 is not defined and the evolution equation 
satisfied by U-{t, x, y) has no left and right limits at t = 0. 
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Case iv) Let us, now, consider the identity 

\y-\Q\"^\ 



X 

From it at t = we have 



and 



5^ 



|12xi|i/2 



2 X 



dx 



G-it,x,y) 



_ 1 sgn y 

t=0 TT X 

1 Sgn?; 



t=o 



By deriving ( ^.60 ) with respect to t we get 



d_ 
di 



8^ 
dy"^ 



6 In 



1/2 1 



\y-\Q\ 

|12a;t|i/2 



and, then, by deriving with respect to x we deduce that G_ satisfies ( 2.21 ). 
We have therefore 



and 



U.{0,x,y) = — 



— OO <J OO 



dy' dx'- 



1 



— OO 'J OO 



(,_,7)lC^o(0,x,y) 
]dy' (dx'—^Uo{Q,x,y). 

oj J X- X' 



Note that the initial data satisfy the constraint 
In aU four cases 



JdxU{t,x,y) =0 fori^O. 



(2.60) 

(2.61) 
(2.62) 

(2.63) 



(2.64) 
(2.65) 

(2.66) 
(2.67) 



In the cases ii) and iv) this condition is satisfied also at t = and, therefore, it can be 
considered as a constraint imposed to the initial data that is conser ved a t all times. On 
the contrary in the cases i) and in) the initial data do not satisfy ( 2.67 ), which can be 
considered as a condition on the solution dynamically generated at all times t 7^ by the 

evolution equatio n. 

By inserting (2.26) into ( 2.22| ) we deduce that r can be rewritten as 



{t,x,y) = d^pe-'P^-+'P-y-''P'^H{t,p)Voip) 



(2.68) 



where H(t,p) is the Fourier transform of G{t,x,y). The computation of H{t,p) allows 
us to state the connection with the general formula (|2.11 ) given in the previous section 
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and, specifically, to state the behaviour of V{p) at pi = and at p2 = in the four cases 
considered via the formula 

V{t,p) = -ipi e-'''P"H{t,p)Vo{p). (2.69) 
It is convenient to compute H{t,p) as a limit according to the formula 

H{t,p)^ lim lim +.x^)e*''i^-'^il^l lim /d?/e"*f=^"'=^l^lG(i, x, y). (2.70) 

In fact the distribution {x'^ + x'^)e~'^^^^^e~'^^^'^^G{t, x,y) converges to G{t,x,y) as /i ^ 0, 
ei ^ 0, £2 — > in the sense of distributions and the Fourier transform operator commute 
with the limit. The insertion of the term x'^ + x^ is unessential in the case G± but in 

the cases G± it takes into account explicitly that we have to compute the principal value 
of the integration over x. 

In all cases we make the change of variables 12xsgn< x and y —\xt\^/'^y. For 
the Bessel functions of different kinds introduced in the following we use the notations 
of @. 

Case i) For G+ we get 

H+{t,p) - lim ^ /""da;e^(?i''snt+..i). r'^j^(l_y2)-i/2g^p.|xt|V=j, ^ 

£1^0 TT Jo J-1 

limsgnt/ dxe'^Ti'^gnt+^^O^jg (p2|a;i|i/2) ^2.71) 



and, therefore. 



Case ii) For G"^ we get 



£1^0 

and, therefore 



H+(t,p) = — i— (2.72) 
Pi + lOt 



dxdye'P'''-'P'y^^^Q-^^'^{t,x,y) = 

Jo J-1 

lim ^\t\'/' r^e^(fi^gnt+,£,).j^(p2|^^|i/2) ^ _^ ( i _ ^-^^4] (2.73) 
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Case iii) For G- we get, analogously, 



1 p 

\Pl\ 



(2.75) 



Case iv) The computation of is mucli more involved and the insertion of the term 



becomes crucial. We obtain 



H- (t, p) — lim lim - 



(2.76) 



where 



Ki (-z(p2 + i0)|a;i|i/2) - (i{p2 - iQ)\xt\^/^ 



-1/2 



(2.77) 



and 



/W(P2,2:,<)= lim / dyfe^(f^+-=^)l-*l'''2'-e-'(f^-"^)l-*l'''^)y(y2 + l)-i/2^ 

-\ [Ki (-{P2 + *0)|xi|i/2) + K, (-(p2 - *0)|xi|i/2j _ 
Ki ({p2 ~ z0)|xi|i/2) - ((p2 + «0)|xi|i/2)] . (2.78) 
We need, then, to sum six integrals of the form 

rdxx^-'/'e-'--K,i2PV^) = E^t±J^)IMa-^^^^^^,{f, + l,2;^) (2.79) 
Jo 2/i a a 

extended to complex values of a and /3. Note that ( 2.79| ) a,t p, — 0, due to the presence 
of the gamma function r(/i), has a simple pole which must cancel out by summing up 
all terms in ( 2.76 ). 

In order to perform correctly the analytical continuation in a and (3 and in order to 
extract the singular term at /i = it is convenient to rewrite the confluent hypergeometric 
function of second kind ^! in terms of the confluent hypergeometric function $ which is 
an entire function by using the formula (see (63) at page 261 in p9|) 



r(^)r(^ + l)^{fi + 1, 2; z) = r(^ + l)$(/i + l, 2; z) logz + 

r(A* + i)E 



^-^^j^ [^(m +l + n)- V^l + n)- VX2 + n)]-+ (2.80) 
(2)„ n\ z 
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Therefore we get, for complex values of a and /3, 

At_L2a~'^!^<i>if, + 1, 2; ^)[21og/3 - log«] 
2p a a 



-a 



2/3 

r(M + i 



E 



n— 



+ 1 + n) - V^(l + n) - ^(2 + n)] 



2/3 



-a 



(2.81) 



By inserting it into the right hand side of (2.76) we verify that the coefficient of T{^) at 
fj, — is zero. Therefore the limit /i — )■ can be computed by using the Hospital rule. It 
is convenient to consider first H^{t,p) for pi ^ 0. By recalling that 



we get (pi ^ 0) 



$(l,2;z)--(e^-l) 

z 



H^(t,p) ^ — sgnpie PI 

P2 



The value of (t, p) at pi = must be computed separately and we get 



(2.82) 



(2.83) 



H-it,p) 



Pi=0 



(2. 



We get, therefore, for r 
Case i) 



T+{t,x,y) = id pe 



l2„ „-ipix+iP2V~'itpl-3it^ * 



pi + iOt 



Voip) 



(2.85) 



Case ii) 



P2 + tfJV 



(2.86) 



Case iii) 



T.{t,x,y) = /d2pe-'''^^+'^^^-'*^'-3''*^^Fo(p) H + 0) 



(2.87) 
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Case iv) 

?-{t,x,y) = [d^pe-'P'^+'P^y-^'P'-^^'i ^li^Voip). (2.88) 

J P2 

We are, now, able to give more specific information on the general distributional formula 

T{t,x,y) = f dx'd^Uit,x',y) (2.89) 

J —too 



derived in the previous section. In the cases ii) and it;), since the condition (2.67) is 
satisfied at all times including i = 0, any other operator can be chosen instead of 
/-too ^^'^ ^ and, consequently, dtU{t, x, y) have a definite limit at t = 0. In the case 
i) T has definite right and left limit aX t — Q and dtU{t, x, y) is in general discontinuous 
at t = 0. In the case Hi) r is not defined at t = and at this point the distribution 
dtU{t,x,y) is not regular. 
For V we get 



Case i) 



Case ii) 



Case iii) 



Case iv) 



V+{t,p)=e-'*P'^-''*^Vo{p) (2.90) 



V^{t,p) = -e-*^'?-^^*^ _JEL^Voip) (2.91) 
P2 + iCrO 



V-{t,p) = -isgnpi e-**Pi-3**^ Vo{p) (2.92) 



V-{t,p) = -ie-"p'~3**^ ^Vo{p). (2.93) 

P2 

By taking into account that Vo{p) is supposed to be continuous at pi = and at p2 = 
the equations above show explicitly the behaviour of V{t,p) at these two special points 
in the four considered cases. 

2.2 More on the case of smooth initial data 

We consider here, in more details, the case of initial data U {x, y) (and then V{p)) belong- 
ing to the test function Schwartz space S in the x and y (correspondingly p) variables, 
i.e. the case i) of the previous section. 
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Let us recall that in this case the solution U{t,x,y) of the linearized KPI equation 
satisfies the evolution equation 

dtU{t,x,y)^-dlU(t,x,y)+?,dl ( dx'U{t,x',v) (2.94) 

J —toe 

and can be written as 

U{t,x,y)^ JJdx'dy' Jit,x~x',y-y')U„it,x',y') (2.95) 

where we have introduced for convenience the notation 

J{t, X, y) = d,G+{t, X, y) = !i^a,(12a:t - y'')-''\ (2.96) 

TT 

The initial value of the evolution equation 

dtUo{t,x,y)^-dlUo{t,x,y) (2.97) 
satisfied by Uo{t,x,y) coincides with the initial value of U{t,x,y), i.e. 

Uo{0,x,y)^U{x,y). (2.98) 
Since, as we proved in the previous section, 

J{t, X, y) = j j d^p exp ^—ixpi + iyp2 — Ht — ^ (2.99) 
the solution can also be rewritten as 

U{t,x,y) =y"d2pe-*''^"+'^^''-^*^'-'**^F(p) (2.100) 

where 

V{p) - JJdxdy e'P^^~'P-yU{x, y). (2.101) 

The last formula for U{t, x, y) can be used to get an integral of motion. In fact by 



integrating (2.100) first in y and then in x we obtain that 

fdxJdyU{t,x,y)^V{0)= JJdxdy U{x,y) (2.102) 
for any t. This result does not contradict the condition 

fdxU{t,x,y) ^0, t^O, (2.103) 
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that we proved, in the previous section, to be satisfied at any time t ^ Q since the order 
of integration in ( 2.102| ) cannot be exchanged. 

By means of ( 2.95 ) and ( 2.96| ) we derive that U{t,x,y), for t {t ^ 0) and y fixed, 
decreases rapidly for x —too while for x — > too it satisfies the following asymptotic 
behavior 



U{t,x,y) 



-1 



47rv3te|x| 
1 



327rk|V3tea;2 



dx'dy' U{x' ,y') (2.104) 
dx'dy' [Utx' + iy- yT] U{x', y') + o{\x\-'/^). 



This asymptotic expansion is differentiable in t and twice differentiable in y and obeys 
(^.1|). Therefore we can differentiate in t the condition (2.103) and use ( ^.94 ). This 
procedure leads to the next condition 



dx xUyy(t,x,y) — 0, t ^ 0. 



(2.105) 



Thanks to the asymptotic behavior of U{t,x,y) at large x derived in ( ^.104 ) one y- 
derivative can be extracted from the integral getting 



dx xUy{t,x,y) — 0, 1 7^ 0, 



(2.106) 



but it is impossible to remove the second derivative since JdxxU{t,x,y) is divergent. 
This procedure can be continued to get an infinite set of dynamically generated condi- 
tions. 

Equation (2.104) explains the role of constraints, i.e. conditions of the type ( ^.103 ) 
and ( ^.106 ) imposed on the initial data U{x,y), in the asymptotic behavior of U{t,x,y) 
at large x. Precisely, for each additional constraint that is satisfied we get an additional 
x"^ factor in the decreasing law. 

We want, now, to prove that the dynamical system described by the evolution equa- 
tion (2.94) is Hamiltonian. 

We define the Poisson brackets of two functionals F and G which are differentiable 
with respect to the smooth initial data U{x,y) as follows 



{F, G} = / d^dr, 



SF d 



SG 



Then 
and 



{U{x, y), U{x', y')} = d,S{x - x') S{y - y') 



{F(p), V{p')} = 



-ipi 



SiPi + P'l) Sip2 + P2) 



(2.107) 

(2.108) 
(2.109) 
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Due to (2.97) and ( 2.9§| ) we easily find that, at equal times, also Uo{t, x, y) are canonical 
since they satisfy the same Poisson brackets 

{[/o(i, X, y), Uo{t, x', y')} = dj{x - x') 5{y - y'). (2.110) 

By considering U{t, x, y) as given by ( ^.100[ ) and by using ( ^.109] ) we have 

{Uii,x,y), Uit,x',y')} = 



^d^p jd'p'e 

'ipi 



a I ~ipix+ip2y~itpl~3it^ -ip\x' +ip' ^y' -itp'\-3it'^ 



5{P1+P'i)5{p2+P'2)^ 



(27r)2 

dx5{x - x') 5{y - y'). 
Let us introduce the Hamiltonian 



(2.111) 



= 2 y j'^'^'^y ^' y))' " i y I ^' I y jdxdy{dyU(t, X, y)) {dyU{t, x-x',y)). 

(2.112) 

Notice that the integration in x' in the second term cannot be performed before the 
integration in x and y since the corresponding integral does not exist. 
We prove the following 

Proposition 2.2 

i) The Hamiltonian H is an integral of motion. 



a) The time evolution equation (2.94) Hamiltonian, i. e. it can be rewritten as follows 

dtU{t,x,y)^{U{t,x,y),H}. (2.113) 



Hi) The Poisson structure introduced in (2.107) is not continuous in t at t = 

lim^{U{t, X, y),H} ^ {U{0, x, y),H} 



(2.114) 



Proof. 



(2.115) 



i) Due to ( ^.lOOD definition ( ^.112[ ) can be rewritten as 

H = 27T^Jd^p {pI + 3pIp^^) Vip) Vi-p), 

which is explicitly time independent. Notice that Pi^ in this expression, which has been 
obtained as the Fourier transform of \x'\, is a distribution defined as follows (see |18 ) 

fdp,p^'^ip,)^^raf dp,^^Ell_m^ (2.116) 



Pi 
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where f{pi) is an arbitrary test function belonging to the Schwartz space S. 
a) We use ( 2.10g| ) and (2.115) to write down 

{V{p),H} = -^^±^V{p), 



(2.117) 



where 1/pi has to be understood in the sense of the principal value. Now from (2.7), by 
noting that due to the presence of the exponential term for t ^ 



lim V{t,pi,p2) = 
pi— >o 



(2.118) 



in the sense of distributions in the variable p2, we get the evolution equation (2.1C) for 

any t ^ 0. 

in) Considering the limit for t — s- ±0 we get by using the lemma A.l in the appendix 

lim {Uit, X, y),H}^ -^^{/(O, x,y) + 3 f dx' 8^(0, x\ y). (2.119) 

./too 

On the other side due to ( |2.8| ) 

{UiO,x,y),H}^{U{x,y),H} = 

Id^pe-'P'^'+'P^y {Vip),H} = 

dlU{0, x,y) + lf dx' sgn{x - x')dlU[0, x' , y), (2.120) 



where in the last line we used (2.117). 



3 Nonlinear case 

3.1 Initial remarks and notations 

The analysis of the linearized version of the KPI equation, as performed in the previous 
section, indicates clearly that it is convenient to study the KPI equation in the space 
p = {pi,P2\ of the Fourier transform of the solutions 

v{t,p) = v{t,pi,p2) - ^ J ldxdye'<'P'-P-y^u{t,x,y). (3.1) 

Correspondingly, also the nonstationary Schrodinger equation 

{-idy + dl-u{x,y))^ = Q, (3.2) 

which is the spectral equation associated to KPI, has to be stated in the Fourier space 
of the variables p — {pi,p2}- 
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If $(a;, y, k) is the Jost solution of ( p.2\j the corresponding Jost solution in the p— space 
can be defined by first shifting <i>(a;, k) and then by taking the Fourier transform 
according to the following formula 



(3.3) 



The specific role played by the spectral parameter k is emphasized by separating it from 
the p variable by a | and a bold font is used to indicate that k belongs to the complex 
plane. 

Then the spectral equation (3.2) takes the form 

[£(p)-2pik]i.(p|k) = j<fp'v{p-p')v{p'\k), (3.4) 

where 

C{p)=P2-pI. (3.5) 

The integral equation determining the Jost solution 1^ has, in the Fourier space, the 
form 



where 



p{p\k)= jd^p'v{p~p')v{p'\k) (3.7) 

and 5^{p) = 5{pi)5{p2)- The new function 

p{p\k) = [C{p) - 2pik] v{p\k) (3.8) 

will play a crucial role in the following. Let us introduce the boundary values of v and 
p on the real axis of the k-plane 

v^{p\k)^ lim v{p\k), p^{p\k)^ lim p{p\k), fc = kju (3.9) 

and write 

.^{p\k) = 6\p) + ^y^l . (3.10) 
C[p) - 2pik T «Opi 

Then the spectral data are most naturally obtained (see ||l^) by computing p{p\k) at the 
special value of k that vanishes the denominator 



— \ ^ Id ri v{n~ r))v^\ V — 



2pi / 



d^p'v{p~p')i.^{p' . (3.11) 



We conclude that the variables p = {pi,P2} can be considered as the spectral variables of 
the spectral data. For their introduction in the general scheme of the resolvent approach 
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see |TJ, H and for analogous variables used in the KPII case see . They are related 
to the more familiar spectral variables {a,/3} (see e.g. ^ 0|) by the formulae 

£i-p) C{p) 

-2pi 2pi (3.12) 

pi = a - (3 , p2 = - • 

We denote the spectral data in the variables {a,P} by r^(a,/3). They are related to the 
previous ones by 

r±(a, /3) = -2.zp±(/3 -a,(3'- a'\a) , p± ^) = f r± ^ 

V 2pi / 27r V 2pi -2pi 

(3.13) 

Spectral data using the standard spectral variables {a, P} evolve in time according to 
the formula 

r±(t,a,/5) =exp(4ii(a3-/33)) r±(a,/3), (3.14) 
while the spectral data evolve in time as follows 

i.e. in the same way as the linearized KPI equation in the Fourier transformed space. In 
fact this is the deep reason of the privileged role played by the {pi,P2} variables in all 
the theory. 

Note that the Jacobian of the transformation {a,/3} {pi,P2} 

dipi,P2) „| ^1 

-^^ = 2\a-(3\ (3.16) 

is not bounded and not always different from zero. Therefore, the use of the two sets of 
spectral variables can be not equivalent. In fact (see |l^]) the use of the {a, /?} variables 
requires to perform a subtraction in some formulae of the direct and inverse problem, 
while, as we will show in the following, this is no more necessary if one uses the {pi,P2} 
variables. 

Let 7V^ denotes the integral operator with kernel 

niia, P) ^ 6{a ~ 13) T OiT<j{a - (3)) r" {a, (3) , cr = +, - . (3.17) 

Then the characterization conditions for the spectral data can be written as (cf. [|l^ ) 

71171^"^=!, 0- = +,-, (3.18) 
ninl^ =7l'L7V'J . (3.19) 
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If the integral equation in (3.6), (3.7) has no homogeneous sohitions, which we always 
suppose, the operator 7?.^ satisfies also the equation 



-o- 1 Tpo- 



711=1, 



(T = +, - , 



(3.20) 



which is complementary to ( 3.18 ). Two other sets of spectral data are used in literature, 
the so called unitary «S-matrix 



s = nl^ 7^: 



and the self-adjoint integral operator 



:f=f = x. 



(3.21) 



(3.22) 



Notice that, thanks to ( 3.19 ), the operator as defined in ( 3.22 ) is independent on the 
signs ± of the constituent operators 7?.^. Multiplying (3. IE) by 'RZ^'^ from the left and 
by 7^I'" from the right we get iV^'RrJ' = 'Rr^"'^'lV_, due to ( ^.20] ). This proves the 
independence of the <S-matrix defined in ( ^.21 ) on the sign a. Properties of S and 
given in (3.21) and (3.22) follow as well from (3.18) and ( |3.20 ). Taking into account the 
triangularity property of the operators 7?.^ the six equations ( 3.18 ) and ( 3.19 ) reduce to 
the four ones [fTl 



/3 

r/3 \ 



r'^{a, 7) r'^ (/3, 7) 



and, analogously, instead of (3.2C) we can write 



r'^(a,/3) +r-'^(/3,a) =a dj r"^ {j, f3) r-'' {'j,a) . 

J a 

Using now ( 3.13| ) we can rewrite ( p.23 ) and (3.24) in terms of p 



(3.23) 
(3.24) 

(3.25) 



-2pi 

2iTTapi d^p' sgn{pip[)i^{~acr'{pi-p[)p[) 6{piC{p') - p[C{p)) x 



P' 



2p[ 



P -P 



C{p'-p) 



2{p[-Pi) 



(3.26) 



Notice that due to the (5-function and (3.5) in the integrand on the r.h.s. we have 

£(y)^Ap) C{P'-P) _ C{-p) 

2p[ 2pi ' 2{p[~pi) -2pi ■ 



(3.27) 
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Correspondingly from (3.25) we have 



2inapi hfp'^{{pi~p\)p\) S{ipi~p[)Cip)-piC{p^p')) x 



P \P 



C{p~p') 



2(P1-Pi) 



(3.28) 



3.2 Properties of Jost solutions and Spectral Data 
3.2.1 1/k expansion of the Jost solution at t — 

To get the asymptotic 1 /k-expansion of the Jost solution we consider the expansion of the 
kernel {C{p) — 2pik)~^ of the integral equation (3.6). By considering it as a distribution 
in the variables L{p) and pi we have that 



C{p) - 2pik 



E 



c{py 



(2k)"+i(pi±io/:(p))"+i ' 



oo, ± kcv > 0. 



(3.29) 



Note that the coefficients of the expansion depend on the half-plane in which k tends to 
infinity. Consequently, also the expansion of the Jost solution at large k depends on the 
half plane in which the limit is performed. If we note the coefficients of the asymptotic 
expansion as follows 



oo ^ / s, oo 



n=0 



(2k)" 



k oo, ± kcj > 



(3.30) 



we have from (3.6) that 

^oip) ^ S'^ip) , Po{p)^'"{p)^ 



<{P) 



E 

m=0 



[pi±iOL{p)Y 



n = 1,2,, 



where Pn{p) due to the (3.7) and (p. 29]) obey the recursion relation 



Pt{p) 



ri-l 

E 

m=0 



(P qv{p~q)C{qy 
[qi ± i{)C{q)Y 



■Pm(<7): 



71= 1,2, 



(3.31) 



(3.32) 



Only the leading coefficients are independent on the sign of kg. The coefficients with 
n = 1 



^tip) 



v{p) 



Piip) 



pi ± iOC(j)) ' 
will be of special use in the following. 



(fq 



v{p-q)v{q) 
qi ± iOCiq) 



(3.33) 
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3.2.2 Properties of Spectral Data 



In order to study the properties of the spectral data we, first, substitute ( p.6[ ) into (3.7) 
getting an integral equation for p(p|k) 



P(p|k) = v{p) , , ^^^^ _ ^^^^ 
Performing the limiting procedure ( [3.9[ ) we have 



P {p\k) = v{p) 



v{p - q) p{q\^) ■ 



vip-q) p'^(qlfc) 



(3.34) 



(3.35) 



C{q) - 2qik T Mqi 

and, then, recalling definition ( 3.11 ), the following representation for the spectral data 

(Pq 



2pi J 



v{p) +pi 



PiC{q) - qiC{p) =F iOpiqi 



v{p-q)p U 



C{p) \ 
2pi ) 



(3.36) 



2pi 



We deduce that for v{p) vanishing rapidly at large p also the spectral data p^yp 
vanish rapidly at large p. However, the spectral data are discontinuous at p = 0. In fact 
if the limit p — )■ is taken along the line p2 ~ 2/3pi, with (3 an arbitrary constant, we get 
from (3.36) 



lim p 

P2 = 2^^P1 



|^)=P^(0|/?)=<0) 



cPqv{q)p^{q\l3) 
C{q) - 2/3gi T ^Oql 

i.e. the limit depends on (3. In addition if we fix p2 7^ we get 



lim p 



^P) \ 
2pi J 



viO,P2) 



and combining the results { ^1% and (^!38|) 



lim lim p 



^r^l = fim lim p^ (pi,2(3pi 

2pi J /J^oopi^O V 



Pi 

2 



viO). 



(3.37) 



(3.38) 



(3.39) 



Note that from ( 3.381 ) it results that the spectral data at pi =0 are smooth (Schwartz) 

r±(/3,/3) = -2^jp±(0|/3) 



functions of P2- Moreover, since 



and 



lim — r 

p^oo 27r 



= lim p 



Cip) \ 
2pi J ' 



(3.40) 



(3.41) 



we deduce from (3.37) and ( 3.381 ) that (no matter how rapidly decreasing at large p is 
chosen v{p)) the spectral data r*(a, /3) do not vanish, in general, at large distances in 
the {a, (3} plane. The spectral data vanish at large distances only if the constraint 



Vi0,p2)^0 



(3.42) 



is satisfied. 
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3.2.3 1/k expansion of the Jost solution for t 

Since the pro perties of the Jost solution depend on the specific form of the spectral 
equation (3^) we would expect that the introduction of a parametric dependence on t in 
v(j)) will not change the asymptotic 1/k-expansion of i^{p|k) or, more precisely, we would 
expect that the only difference would be a parametric dependence on t of the coefficients 
of the expansion. However, we will see that the time evolution of a solution v{t,p) of the 
KPI equation has special properties in the neighborhoods of t = and we have to take 
into account that not v{t,p) but v{t,p) as defined in the following equation 



v{t,p) — e ^ PI v{t,p) 



(3.43) 



is continuous at p = for t ^ 0. 

Therefore in writing the integral equation (3.6) for the Jost solution it is convenient 
to factorize explicitly the singular behaviour in time as follows 



u{t,p\k)^S\p) 



1 



We get 



C{p) - 2pik 
lim k[i^(t,p|k) — (5^(p)] =v{t,p) lim 



(Pqe «i v(t,q) iy{t,p — q\]i.) . 



C{p) - 2pik 



and using proposition A.l in the appendix 

-1 



lim ]i[v{t,p\k) -S^ 



2(pi + iOt) 



''v{t,p) 



-1 v{t,p) 



Pi 



iOt 



(3.44) 



(3.45) 



(3.46) 



Comparing this result with (3.33) we see that the coefficients of 1/k in the asymptotic 
expansions coincide only if 

sgnkQ = sgnitCip)) , (3.47) 

or in other words only under this condition the limits k ^ cxd and t ^ commute. 
Notice that this condition is formulated in terms of the variables p corresponding to the 
fact that the special behaviour of solutions and Jost solutions in the neighborhoods of 
t — are manageable only in the transformed Fourier space. 



3.3 Inverse problem and time evolution 
3.3.1 Formulation of the inverse problem 



By working in the general framework of the resolvent approach |14, 8| we obtained in [|_5| 



the following coupled system of two linear integral equations for the boundary values of 
the Jost solutions on the real fc-axis 



,^^{p\k)^S\p) + 



2 J C{-q) + 2kqi±iqiQ 



E 

r=+,- 



-2^) " 1^ + '^ 



2qi ) 



(3. 
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This formula is the main tool to be used for solving the inverse problem, i.e. for recon- 
structing the potential v{p) starting from the spectral data and for proving that spectral 
data evolving in time as indicated in ( 3.15| ) generate potentials v{t,p) that are solutions 
of the KPI equation. 

In this respect two properties of this system are crucial. First, the system is closed 
in the sense that the spectral data are reductions of the Jost solutions themselves 
via equations (|3.8|), ( |3.9| ) and ( 3.11 ). Second, it depends linearly on the spectral data, 
in contrast with the usual equations written in the literature that are quadratic in the 
spectral data. To our knowledge the only alternative integral equations for solving the 
inverse problem, which are linear in the spectral data, are given in [ p2[ , but they are 
not closed in the sense that they are written in terms of the so called advanced/retarded 
solutions. 



We suppose that the system ( 3.48[) has unique solution in the considered class of 

off~ 



spectral data. Then the r.h.s 
k-plane as follows 



can be continued analytically in the complex 



iy{p\k)=S'{p) + - 



2 J C{-q) + 2kqi 



2qi J 



(3.49) 



For getting from (3.49) an integral equation for p(p, k) it is convenient to use equation 
(B^) computed on the real axis and to rewrite the Jost solution in the r.h.s. of (3.49) as 



ly \^p + q 
Then it is easy to check that 

-'{p + i\t^) _ 1 



2qi ' qiC-iv) +Vi^{-q) - icrOqi{qi +pi) 



Ci-q) + 2kqi C{p) - 2kpi 



7." + £(-.) + 2k,, 



(3.50) 



(3.51) 



where the ratio pi/gi multiplying v in the r.h.s. is well defined due to (3.50). By inserting 
this expression in (3.49) and by recalling definition (3.8) we get the following integral 
equation for /o(p|k) 



p(p\k) 



2qi ) 



2 J C{-q) + 2k(7i 



cr=± 









2gJ^ 




2qi J 



(3.52) 



It has a kernel similar to the kernel of the integral equation ( p.6[ ), (3.7) obtained for J^(p|k) 
in the direct problem. Therefore its asymptotic expansion at large k can be computed 
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by using the same methods used in the previous section. In particular we get that the 
first term in the r.h.s. due to ( 3.30[ ) and ( |3.33 ) is just the potential 



Pi 
2qi 



cr=± 



(3.53) 



that is therefore explicitly reconstructed in terms of the spectral data and Jost solutions. 

It is instructive to consider in details the properties of the function v{p) defined by 
(mi) in the neighborhoods of p = 0. This analysis throws some light on the meaning 
of the characterization equations for spectral data reported in section 3.2. In fact, if we 
rewrite (3.53) by using ( |3.50D as 



v{p) 



1 



-2pi 



2gi 



^ qiC{p) +piC{-q) - iaOqiiqi + pi) 



(3.54) 



and perform the same limit as in (3.37) we get 

1 f- 



lini -D(pi,2/3pi) 

pi->±0 



(fq 



-(0|/3)- 
C{~q) + 2(3qi T iaO 



(3.55) 



This result seems to contradict the fact that the potential v{p) was chosen to be contin- 
uous at p = 0. However, the term in the r.h.s which depends on the sign of pi is y'^'^o- 
portionaltoX;CT=±o-/c^^g'5(£(-q) + 2/3gi)(/F°'p'^)(^g|^^ and, then, due to (|1|) and 

O to E.=± ^Sdar (a, /3)r- (a, /3) = E.=± ^ ( {^^7 ^ - 7^1" ^ ) (7?:^ - ) ) (/?, /?) . 

But due to (|^ and (|^ ("^T"^ - nZ''^){7l\ - TVL) = 22 - S - which is cr- 
independent. Therefore this term is equal to zero and the limit does not depend on the 
sign of pi . The independence on f3 can be proved by using ( [3.37 ) . 

From (3.49) by using the expansion (3.29) we can get the coefficients of the asymptotic 
expansion at large k of J^(p|k). Of special interest is the coefficients of 1/k. Recalling 
( [S.35 ) we have 



v{p) 



Pi ± iOCip) 



<Pq 



2{qiTiOC{-q)) 



2qi ) 



2qi ) 



For the all set of coefficients of the asymptotic expansion of p(p|k) we get 

d?q{-C{-q)y 



PniP) 



2(<7i T zO£(-(z))"+i 



n-p 



P \P 



CW)\ 

2qi J 



(3.56) 



(3.57) 
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3.3.2 Jost solutions for i ^ 

We choose spectral data to evolve as in ( 3.15| ) 



UW - 



exp —it 



pt+M 
pi 



(3.58) 



where P^{p\^^) denotes the initial value of spectral data determined in terms of v{p) 



and Jost solutions i/^(p|fc) according to formula (3.fl). 

Then the Jost solution at any time is reconstructed just by inserting these spectral 
data into ( 3.49| ) 



^{t,p\k) = 5^{p) 



i9 it - 

d'^qe 



27 C{-q) 



CT — it 



2gi 



(3.59) 



In the limit k — > fc ± zO (cf. ( 3.48 )) we obtain the integral equations for the boundary 
values on the real axis, i.e. for the functions v^{t^p\k). In what follows we always suppose 
that these equations are uniquely solvable. 

We have, first, to pr ove that this function v{t, p\k) is indeed the Fourier transformed 
of the Jost solution of ( |l.2| ) , i.e. that it satisfies ( |3.4| ), and, second, we have to find its 
time evolution equation, which is needed to know in order to show that the reconstructed 
v{t,p) satisfies KPI. 

3.3.3 The recontructed spectral equation 



First we prove that z/ defined by ( 3.5£ ) obeys the spectral equation (3^) at any time. 
Let us define p{t,p\\i) by (3.8) rewritten for any time 

p{t,p\k) = [Cip) - 2pik] iyit,p\k). (3.60) 

Then, by the same procedure used to get (3.52) from ( |3.49| ) we derive from (3.59) that 



p{t,p\k) = v{t,p) + 

where we have introduced 
v{t,p) 



1 f <Pqe' 



2 J C{-q) + 2kqi 



;E^:i:^(C)''"('--'l|r) 



f , Pi ,t£l±!2i ^— / C{q)\ 



V {t,p + q 



Ciq) 



2qi 



(3.61) 



(3.62) 



Now due to the assumption on unique solvability of (3.59) we have from (3.61) that 
p{t,p\k) — f cPqv jt, p — q):^{t,q\'k), i.e. (3.7) is satisfied for any t. Thanks to ( [3.8D this 
proves that (|3.4|) is satisfied at any time. 
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To prove that v defined by ( |3.59 ) obeys the integral equation ( |3.6| ) at any time we 
use ( ^^ , (HH) and ( |3J^ ) to get 



i.(i,p|k)-<52(p) 



1 /■ d qe 



2„ ii- 



2 J C{-q) + 2kqi ^ 



E 



<T = ± 



2gi / 



v^(t,p + q 



^q).£{p + q)-2{p, + q,)^ p-{t,p + q\^ 



2qi J 



C{p) - 2kpi 



(3.63) 



Then ( |3.6| ) at any i follows by using definition ( |3.60 ) of p'^ in the r.h.s. This proves that 
v{t,p\k) in (3.59) is the Fourier transform of the Jost solution of the equation (1.2) with 
potential the Fourier transform of v{t,p) defined in ( 3.62| ). 

To prove that v{t,p) satisfies equation (3.43) with v{t,p) continuous at p = it is 
convenient to insert (3.50) (that we proved to be satisfied at any time) into ( [3.62 ) getting 



-,;.£l±M 1 



<T = ± 

'(fqpie "1 2_\ 



C{~P) 



-2pi 
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cr=± 



qiC{p) +piC{-q) - iaOqi{qi + pi) 



(3.64) 



The time evolution of p"^ ^t,p 



in the r.h.s. is not explicitly known. However, it 

is useful t o fac torizc the time evolution it would have in the limit v{p) by writing 
(see also (3.81) below) 



pit,p\k) = e-2«*[3(k+Pi)£(rt+2p?] p(^t,p\k) . (3.65) 
Note that for the special reduction furnishing the spectral data we have that 



P^{t,q 



2qi J 



is time independent. 

Finally we can write by using the identity 



2qi J 



(3.66) 



'1' + ^'^"^ -.6jC:{p)^-2[3p,C{p) + 2pl] 



2qi 



+ qi+pi\ L{p + q) + 2{pi+qif 



(3.67) 
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that 



v{t,p) - v{p) 



qiC{p) +piC{-q) - iaOqi{qi + pi) 



(3.68) 



exp 6itC{p) 



2pi 



2gi 



P''[t,p + q 



At the special values of p and q for which the denominator in the r.h.s. vanishes 

p + q _ _ = {p-a,/3^ -a^] , (3.69) 



with 



gi£(p)+pi£(-g)=0 



a 



^ , „ (3.70) 

and we obtain that, due to (|3.13|) and (|3.66|) , p and p reduce to the same spectral data 

p'"(/3-a,/32 -a^|a). (3.71) 

Therefore, at qiC{p) + piC{—q) ~ the function in the square brackets in the r.h.s. 
vanishes and pi can be extracted from the integral. We conclude that 



lim v(t,p) = u(0). 

p—*0 



(3.72) 



The continuity of v{t,p) at p = ensures us that the exponent in ( 3.43| ) is the dominant 
term in the time evolution of v(t,p) and therefore we recover for u{t, x, y) at large x the 
behaviour computed in the linearized case in (2.104). 

3.3.4 1/k expansion in the inverse problem 

One can compute the coefficient of 1/k in the asymptotic expansion at large k of i^(i,p|k) 
also by using the Jost solution as reconstructed from spectral data via the integral equa- 
tion ( 13.591 ). 

Thanks to proposition A.l in the appendix we have that 



lim k[iy{t,p\k)~S^{p)] 



1 fd qe "1 \ 



qi — iOt 



cr=± 



As mentioned in the appendix the iOt term is not relevant for t ^ and, therefore, we 
get from (3.62) 



lim k[Ki,p|k)-<52(p)] 

k^oo 2pi 



for t 7^ , 



(3.74) 
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The singular function 1 /pi does not need any regularization at t 7^ since it is smoothed 
by the oscillating behaviour in time of v{t,p) indicated in (3.43). On the other side at 
i = (3.73) is discontinuous 



lim lim k\i^{t,p\k) ~ S^{p)] = 

t^cOk^oo 



<7 = ± 



(fq — 



qi — ieO 



p-" q 



2gi ) 



p + q 



C{q)_\ 
2<Zl J 



(3.75) 



Substituting (gi — ieO) ^ = {qi T *'C(— (7)0) ^ + 2ine5{qi)d{±€q2) we can use (3.56) and 
obtain 



v{p) 



lim lim k[iyit,p\k)~SHp)]=- 



~2~ 



<T = ± 



lim P- 



2gi y 



2(?i 



^?(±e(72), e = +,-.(3.76) 



Due to ( 3.38 ) and ( [3.50 ) the last term is equal to iTre'&{^ep2)S{pi)v{0, P2) and, conse- 
quently, 

lim lim k[Kt,p|k) - S\p)] = -T^r^^^^ , (3.77) 



2{pi + ieO) 



Comb ining this result with the previous formula ( 3.74 ) we recover the same coefficient 
( ^.46 ) computed by usi ng th e integral equation for i/(t, p|k) of the direct problem and, 
by recalling ( 3.73 ) and ( 3.74 ), we obtain that 



Pi + m 



1 fd^qe 



i2 it- 



qi — iOt 



-i:^(«lt^)^1'..+.|t^) 



(3.78) 



3.3.5 The recontructed time evolution of the Jest solution 

To derive the time evolution equation of the Jost solution we differentiate with respect to 



t equation ( 3.59 ). Proposition 
of this equation 



in the appendix furnishes the derivative of the kernel 



dZT— — 
e "1 

dt C{-q) + 2kqi 



= 3i£(-9)- 



qi — iOt 



2i{2qf + 3iqi-k)Ci-q)) 



By using ( |3.78| ) and the identity 



^-^=6Cip)(k+^ 
91 V 2(71 



C{-q) + 2kqi 
-2[3{k + pi)Cip) + 2pl] + 



(3.79) 
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'{^ + 91 + PI ) + + 2(Pi + 



(3.80) 



we obtain an integral equation of the type (3.59) for the quantity di'{t, p\\i) / dt+2i 1^(1, p\k) 
[3(k+j3i) C{p) + 2pf\ with the inhomogeneous term equal to —3iv{t,p) -^^^^^ Thanks to 
the assumption of uniqueness of the solution of this equation we get, finally, the searched 
time evolution equation for the Jost solution 



dt 



= -2i:/(i,p|k)[3(k + pi)/:(p) +2p3] Jd^qv{t,q) 



-^-(^,P"9|k). 
qi + zOr 



As expected the time derivative of i^(p|k) is discontinuous at i = 



(3.81) 



lim 



di^it,p\k) 
dt 



-2il^{p\k)[3{k + pi)C{p) + 2pl] ~3l (fqv{q) 



qi±iO 



v{p - g|k) 



(3.82) 

Note that at large k the leading term in the l.h.s. is of order 1/k while in the r.h.s. it is 
of order k*^ and, due to ( 3.33 ), it is equal to 



m_ 

pi + iaC{p)Q pi ± iO 



a — sgn ko 



(3.83) 



This term, if the initial data v(j)) do not satisfy the constraint v{0,p2) = 0, is zero only 
if sgnkcj = sgn(i£(p)). Once more we recover the relevant fact that the two limits t 
and k — > cxD can be exchanged only if condition (3.47) is satisfied. 

3.3.6 The reconstructed KPI equation and the integrals of motion 



Recalling (|3.8D, that we proved to be satisfied at all times, we deduce from (|3.8lD that 

■ v{t, p) - 2zpi [iC{p) + 2p2] p{t, p\k) - 
3»[£(p) - 2pik] [d\v{t,q)-^^ [v{t,p-q\k)-5\p-q)\ . (3.84) 



2,1 ■ 



dt 

"pi + iOt 



The time evolution equation of v{t,p) can be obtained from this equation by computing 
the coefficient of the 1/k ter m in its asympto tic e xpansion at large k. 

We first note that from ( 3.46| ), by using (3/7) that we proved to be satisfied at all 
times, we get 



lim k[p{t,p\k) - v{t,p) 



2 v{t,p- q)v{t,q) 



d^q 



2(gi + iOt) 



(3.85) 
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Then we multiply (3.84) by k and compute the limit for k — > oo by inserting ( |3.46 ) 
and ( ^!85| ). We obtain 



dv{t,p) ^ 
dt 



\Pi + iOt 



iOt 



+ 6piC{p)+ipi v{t,p) 



3ipi Id qv{t,q) 



C{q) v{t,p-q) 



qi + iOt pi — qi + iOt 



(3.86) 



which can be easily transformed into the evolution form of the KPI equation (1.15) 
reported in the introduction. 

Of special interest is the behaviour of equation (3.84) in a neighborhood of the point 
of discontinuity p — 0. Precisely we consider the limit p ^ along a path such that 
limp^oPi/Pi = 0- Note that, because of ( 3.43| ) and the properties of v{t,p) at p = 0, 
only in this limit v{t,0) = v{0). Then all the terms in the r.h.s. go to zero and we get 
that 

hm M!« = o. (3.87) 

Pj/Pl^O 

From ( ^ieil) , recalling (jsissl) , we get the dispersion relation 



C(k) = lim p(t, Oik) = w(0) 

p — >0 

pS/pi — 



(i^q 



2 J C{-q) + 2kq, 









>C(g)\ 




2qi . 




2qi ) 



(3.88) 

We conclude that C(k) generate an infinite number of conserved quantities, that can be 
obtained by the recursion relations (3.32) computed at p = 0, when the limit is taken in 
the specified way. These conserved quantities are expressed in terms of the spectral data 
by means of equations (3.57) evaluated at p = 0. 



4 Concluding Remarks 

We conclude by attracting the attention of the reader to some specific features of the 
formulation of the direct and inverse problems considered in this article and in [l^ . 

1. The direct and inverse problems are formulated in terms of the Fourier transformed 
Jost solutions. It results that the corresponding variables pi and p2 are more proper 
for the study of the properties of these solutions than the original ones x and y. 
For instance condition ( |3.47 ) which play a relevant role in the theory cannot be 



reformulated in terms of the x, y variables. 

2. We use variables pi and p2 and their combinatio ns as argum ents o f the spectral data 
instead of the standard variables a and (3 (see ( |3.12 ) and ( 3.13| )). These variables 
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are suggested by the form of the direct problem. Properties of the spectral data 
are more naturally formulated in their terms. 

The main tool of the inverse problem is a set of two integral equations for the 
boundary values of the Jost solutions on the real axis of the complex plane of the 
spectral parameter k. This set depends linearly and not quadratically, as in the 
standard formulation, on the spectral data. In addition the spectral data themselves 



are given in terms of the Jost solutions by means of ( 3.11 ) . This closed formulation 
linear in the spectral data enables us, for instance, to study easier the properties 
of the Jost solutions under time evolution and t o deduce easier the asymptotic 
expansion at large k and dispersion relation ( 3.8§| ). 



Functional character of spectral data and their different properties according to the choice 
of spectral variables require some special remarks. 



If we choose the standard set of spectral data ( 3.13 ) we get spect ral d at a whi ch have 



bad behaviour for a and (3 going to infinity. More precisely, due to ( ^3.37 ), ( 3.38 ), ( |3.40 ) 



and ( 3. 41] ), the diagonal values r^(/3, 13) and values 'close' to the diagonal do not decrease 



at infinity (see also fl^]). It is necessary to emphasize that this happens for initial data 
of KPI as small in norm and smooth as we like. The only p ossib ility to remove this 
constant behavior is to impose the constraint ?;(0,p2) = (i.e. ( |l.5| )) on initial data. In 
this case we get a, 1/ (3 behavior, to improve which we need a second constraint and so 
on. Due to ( ^3.37] ) it is clear that in order to get spectral data r='=(/3,/3) decreasing for 
/3 — !■ oo quicker then any negative power of (3 it is necessary to impose an infinite set 
of constraints on initial data. Let us mention that this situation is specific of the 2 + 1 
dimensional case and has no analogue in 1 + 1. 

In ||l^ it was shown that the study of the properties of the Jost solutions by means 
of the inverse problem in terms of spectral data r^(a, [3) needs a special care due to this 
asymptotic behavior. In order to compute the asymptotic 1/k-expansion or the time 
derivative of the Jost solutions it is necessary to use a special subtraction procedure to 
avoid divergent expressions. 

The use of the spectral data {p, ^^^^ article enables us to remove this 

complications. This functions rapidly decrease for p — > cxd in all directions (roughly 



speaking as initial data v{p)) and are smooth at pi ^ (see ( ^.36 )). Of course difficulties 



cannot disappear thanks to a convenient change of variables and in fact also the spectral 
data [p^ belong to the Schwartz space since they are discontinuous at 

p = 0. This difficulty is, however, more naturally manageable and direct and inverse 
problems can be formulated practically in the standard way. Special attention must be 
paid to the discontinuity of the time derivatives at f = 0. This problem is controlled by 
the lemma and propositions given in the appendix. 

As shown in section 3 the properties of all the quantities under consideration are 
essentially different for t — Q and t ^ Q. This situation is specific of the multi-dimensional 
integrable evolution equations, which are non local already in their linearized version. 
In fact the nonlocality implies that the dispersion relation of spectral data is singular as 
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indicated for instance in (3.58). Then for t ^ the rapid osciUations of the exponent 
in ( 3.58| ) near pi = can control a singularity at pi =0 while for t = a specific 



regularization is needed. This is hidden if we use the standard spectral data r (a, In 



fact the exponent in (3.14) changes the functional class of the initial data. The diagonal 
values r^(t,f3,f3) are time independent and thus have the same bad behaviour at large 
/3 as the spectral data at time t = 0, but asymptotic expansions 'close' to the diagonal 



of the type (3.41) show again rapid oscillations for t and p2 different from zero. 
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A Appendix 

Lemma A.l The function exp{iT / pi) defines a distribution in the Schwartz space of the 
variable pi depending continuously on the parameter r. 

i) This distribution is continuously differentiable in r for any r ^ 0. Precisely, we 
have 

5rexp( — ) = — expf — ), r ^ 0, (A.l) 

\PiJ Pi \PiJ 

where the r.h.s. is a well defined distribution in the same space, 
a) At T = there exist right/left limits 

lim — exp f — I = — . (A.2) 

Thus dr exp can be considered a distribution also in the variable t and we can write 

for arbitrary r 

exp ( — ) = ^— — exp ( — ) . (A.3) 

\PiJ pi-iUT \Pi J 

Proof 

i) Let us consider 

/(t)= Jdp, e'^/fV(pi) 
where f{pi) G S. Subtracting and adding terms we can write 

/(r)= /dpie-/Pi[(p(pi)-z9(l-|pi|)^(0)] +(^(0) f dp^e'^'^K 
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Changing pi to 1/pi in the second term we have 

/(r)= /dpie*-/^^[(^(pi)-^(l-|pi|V(0)] +^(0)/ ^e'-^\ 



\pi\>i Pi 

Now it is obvious that both terms are difFerentiable in t for t 7^ and 

drfir) = * e^^/P^ [^(pi) - ^{1 ~ \pi\MO)] + MO) f — e^"'^- 
J Pi J\pi\>iP^ 

Substituting again pi for 1 /pi in the second term we obtain that it cancels out with the 
second term in brackets. This proves i). Let us remark that the last formula shows that 
for r ^ it is not necessary to regularize the factor 1 /pi in front of the exponent in ( A. 1 ) . 
As well one can consider this factor indifferently as principal value or as (pi ± iO) . 
ii) Again subtracting and adding terms we can write 

/— e^^/P^^bi)= / ^e-M [^(pi)_i?(l_|pi|V(0)] + 
J Pi J Pi 

sgnT^(O) / ^e'P\ 
J\pi\>\t\ Pi 

where now in the second term pi was substituted for rpi . To compute the limit for t 
notice that 



lim 

T->0 



dpi 
e ' 

pi|>|t| Pi 



Thus 



lim / ^ e-/^Vbi) = / ^ Mpi) - ^{1 - \pi\MO)] ± ^MO) ^ 
■^^oj Pi J Pi 

— ip{pi)±i'K(p{0), 
Pi 

where to get the second equality we chose p^^ in the integral as the principal value. Then 
( |A.3[ ) follows due to ( A.l) and the remark made after the proof of i). 
Some useful propositions follow from this lemma. 

Proposition A.l The function p^_2kpi defines a distribution in the Schwartz space of 
the variables pi, p2 depending continuously on the parameter t. 

i) This distribution is continuously differentiable in t for any t ^ Q and has well 
defined right/left limits at t = according to the following formula 

d e~^'*5^ 6ikp2 -sit^ 3ip2 sit^ , . 

— -e PI ■ — e PI (A. 4) 



dt p2 — 2kpi p2 — 2kpi pi + iOt 
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lim — ( 

k-too p2 — 2kpi 



.3it^ 1 e "1 



PI = — 



Note that, due to (3.29), 



lim 

k^oo p2 — 2kpi 



1 



2{pi + i0ap2) ' 
and therefore we get from the previous proposition: 
Proposition A. 2 The limits in the following formula 

If 

lim lim e~ 

t^±Ok^oo p2 — 2kpi 

can be exchanged if and only if 

sgnkg. = sgntp2. 



2pi + iOt 



a — sgnkc 



3it- 



From the point m) of lemma A.l we get that 
Proposition A. 3 For any distribution defined by 

Uit,x,y,.Jd^pe~e.,[ 



Pi 



Vip) 



with V{p) an arbitrary function belonging to the Schwartz space 



lim 



lim 



dx'U (i, x' ,y) — / dx'U (0, x' , y) 



=poo 



dx'U {t, x' , y) = dx'U{0, x', y). 

o J =pc30 



(A.5) 



(A.6) 



(A.7) 



(A.; 



(A.9) 

(A.IO) 
(All) 



Therefore the operator inverse of dx defined by J^^^dx' commute with the limit t 



lim 



dx'U{t,x',y) = dx'U{0,x',y) 



=poo 



while the antisymmetrical inverse operator satisfies 

lim — 

t^±o 2 



+ / ]dx'U{t,x',y) ^1 dx'U{0,x',y). 



(A12) 



(A13) 
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